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The coupling between a primordial magnetic field and the cosmic axion field generates a helical 
component of the magnetic field around the time in which the axion starts to oscillate. If the energy 
density of the seed magnetic field is comparable to the energy density of the universe at that time, 
then the resulting magnetic helicity is about \IIb\ — (f 0~^''G)^ kpc and remains constant after its 
generation. As a corollary, we find that the standard properties of the oscillating axion remain 
unchanged even in the presence of very strong magnetic fields. 
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I. INTRODUCTION 

The observations of galaxies and galaxy clusters un- 
ambiguously show that we live in a magnetized uni- 
verse. Large-scale magnetic fields with coherence length 
^ ^ Mpc and intensities B ^ /iG have been observed 
in all gravitationally bound large-scale structures. How- 
ever, whether these magnetic fields have a primordial ori- 
gin or were generated during the period of galaxy forma- 
tion is still an open question (for a full discussion see 
Refs. 1^, iJ, H). Understanding this point is a ma- 
jor goal of modern astrophysics and cosmology, and for 
this purpose it is extremely important to consider all the 
phenomena which leave a different signature on magnetic 
fields of different origin. 

As we are going to show below, one of these phenom- 
ena is the interaction of the magnetic field with a cosmic 
axion field. We will show in particular that the primor- 
dial axion oscillations generate a helical component of the 
magnetic field. Of course, this effect regards only mag- 
netic fields of primordial origin, since it takes place in 
the very short interval of time when axions start to oscil- 
late. This happens when the temperature of the universe 
is about I GeV, much before the galaxy formation. In 
other words, if we believe in the existence of axions then 
any primordial magnetic field observed today must be (at 
least partially) helical. 

The magnetic helicity is a very peculiar quantity asso- 
ciated with a magnetic field Q, and speculations about 
the generation ofprimordial helical niagnetic fields exist 
in the literature In general, 

given an electromagnetic field A'^ = {A^,A) we define 
the magnetic helicity as 



HB{t) = / (fx A- V X A, 



(1) 



which can be considered as a Chern-Simons term be- 
cause of the relation IJ^^ d^xF^^Ft"" = HB{t2)- Heiti), 
where i^^^ = d^Ai, — d^A^ is the electromagnetic field 
strength tensor, F^"' = (1/2^^) e'"''"^-Fpa its dual, 
g = det being the determinant of the metric tensor 



and e^^'^P'' the Levi-Civita (pseudo-)tensor. The defini- 
tion of the helicity given above is valid in the Minkowski 
as well as in the expanding Freedman-Robertson- Walker 
universe, since it is easy to see that Eq. Q is invariant 
under general coordinate transformations. On the other 
hand, it is odd under discrete P and CP transforma- 
tions (it is a pseudo-scalar quantity). This leads to the 
well known result that the presence of magnetic helicity 
in our universe is a manifestation of a macroscopic P and 
CP violation. 

Of course, the phenomenon we are investigating cru- 
cially depends on the existence of the axion field. This is 
related to the so-called Peccei-Quinn (PQ) mechanism 
[l^j l for the explanation of the smallness (or absence) 
of the CP— violating part of the QCD Lagrangian. Af- 
ter almost 30 years, this is still the most trusted argu- 
ment for the solution of the above problem, known as 
the strong— CP problem (for a general review see, e.g., 
Ref. 01). This explains why the existence of the axion is 
largely believed, even without any experimental evidence 
for it, and justifies the great interest in its phenomenol- 
ogy. 

Axions are chargeless, weakly interacting, pseudo- 
scalar particles which emerge as (pseudo-)Goldstone 
modes of the (almost) conserved PQ symmetry J7(1)pq. 
Although chargeless, axions interact with photons by 
means of the anomalous term ga-yuFF, where ga-y = 
Q:em/(27r/a), IS the electromaguetic fine structure 

constant, and fa is the scale at which J7(l)pQ is spon- 
taneously broken, known as the axion or PQ constant. 
This scale characterizes all the axion properties on the 
phenomenological ground. Even though its value is not 
predicted by the PQ mechanism, and thus is model de- 
pendent, it is constrained to the very narrow allowed win- 
dow 10^° ^ fa ^ 10^^ GeV by astrophysical and cosmo- 
logical observations and considerations 0, 0| . This is 
a remarkable point since it will make all our results es- 
sentially independent on the specific axion model. 

It is quite interesting (and even rather surprising) that 
the axion field, introduced to solve the strong— CP prob- 
lem in the QCD Lagrangian, generates a helical compo- 
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nent of the primordial magnetic field, contributing to the 
macroscopic CP— violation in the universe. This point 
can be more easily understood considering the cosmologi- 
cal evolution of the axion field. Axions are born when the 
universe is at the temperature T ~ fa, and C/(1)pq spon- 
taneously breaks. After the universe has cooled down 
enough for the axion mass to be comparable with the 
Hubble expansion rate H , axions start to oscillate driving 
the QCD Lagrangian to its CP— conserving minimum. 
During that time axions constitute a (quasi-)zero mo- 
mentum condensate, weakly interacting with the exter- 
nal magnetic field. The main effect of this external field 
is to "extract" axions from the condensate. This process 
is governed by the simple Eq. 1)28(1 , where Ua is the num- 
ber of axions per co-moving volume and e is a positive, 
time dependent function, proportional to the magnetic 
energy, and inversely proportional to the conductivity of 
the primordial plasma. For vanishing magnetic fields, the 
number of axion per co-moving volume remains constant 
during the time of coherent oscillations. This is a well 
known result (see, e.g., Ref. and references therein) 
which leads to the upper limit on the PQ constant (the 
so-called cosmological bound [l^l). We will see later that 
this bond remains essentially unchanged even in the very 
strong field limit, owing to the smallness of e. 

However, the effect of the axion primordial oscillations 
on the magnetic field is more interesting. In fact, the 
evaporation of axions from the condensate is followed by 
the productions of photons and a resulting increment of 
the magnetic field strength. Since the probabilities of cre- 
ating left-handed and right-handed photons are different 
(CP is not exact), the final result is an overproduction of 
photons of one kind with respect to the other. In other 
words, the helicity of the system changes. 

More specifically, the axion coherent oscillations play 
the role of an a^— dynamo on the evolution of the mag- 
netic field because of the term adynV x B [see Eq. 
where adyn is proportional to e. This term is P and CP 
odd, and is responsible for the generation of the magnetic 
helicity. As we expect from the above discussion, for this 
case the dynamo is not efficient (ctdyn is small), owing 
mainly to the smallness of the axion-photon coupling ga-y 
and to the large value of the primordial conductivity a. 
This last point refiects the known result that helicity is 
conserved when the conductivity is infinitively large. 

The interaction between pseudo-scalar and electromag- 
netic fields has been the object of various papers in the 
literature. 

In the seminal paper by Turner and Widrow j^, it was 
suggested that, during the inflationary epoch, small fluc- 
tuating magnetic fields could have been amplified due to 
the coupling with the axion field. This possibility was 
exhaustively studied in the subsequentpapers by Gar- 
retson et al. M, and Field and Carroll 
Ahonen et al. [2l| studied the evolution of the cosmic ax- 
ion field in the background of an external homogeneous 



magnetic field just after the QCD phase transition, while 
in a recent paper, Lee et al. 22] investigated the possibil- 
ity of generating magnetic fields through the coupling of 
an evolving pseudo-scalar field during the period of the 
Large Scale Structures formation. 

In the next Section we will study more accurately the 
mechanism of helicity generation by primordial axion os- 
cillations before the QCD phase transition. 

II. COSMIC AXION OSCILLATIONS AND 
MAGNETIC HELICITY 

We start with the following interaction Lagrangian be- 
tween the cosmic axion field (j) the electromagnetic 
field Af^ in curved space-time 

+ \ga.,ci,F^,F^'' +fA^), (2) 

where the external current takes into account the in- 
teraction between the cosmic plasma and the primordial 
magnetic field (see below). 

Physics becomes more clear if we introduce the electric 
and magnetic fields E and B. In a flat universe described 
by a Robertson- Walker metric, ds^ = dt^ — a^dx^, where 
a(t) is the expansion parameter normalized so that at 
the present time to, a(to) = 1; this operation can be 
performed in the usual way as 

Foi = —aEi, Fij — eijkO^Bk, (3) 

where Latin indices range from 1 to 3. We shall work in 
the Coulomb gauge A" — diA"^ = 0, in which Eq. J^Jl be- 
come aE = —A and a^B = V x A, where a dot indicates 
the derivative with respect to the cosmic time t, and the 
spatial derivatives are taken with respect to co-moving 
coordinates. 

In terms of the electric field, the external current has 
the form = (0,(tE), where the conductivity of pri- 
mordial plasma is a temperature-dependent function that 
can be expressed as cr(T) = kT. For Aqcd ^ ^ ^ ™w 
(where Aqcd — 200 MeV and mw — SOGeV), k is a 
slowly increasing function of temperature of order unity 
[13 (k ~ 0.76 for T ~ Aqcd, k ~ 6.7 for T ~ mw)- For 
very high intensities of the magnetic field, that is above 
the critical value Be = ml/e ~ 4.4 x 10^'^G {uie and e 
are the mass and absolute value of the electric charge of 
the electron), the expression for the conductivity needs 
to be multiphed by B/B^ |2i |. 

Introducing the "angle" Q = (j)/ fa, the Hubble param- 
eter H = a/a, and neglecting any spatial variation of 0, 
from Lagrangian we get the equations of motion 

e + 3He + m2e- ^E-B, (4) 

. . . 

=J+JLi+Je: (5) 
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where j = crE is the Ohmic current, = E + 2HE is 
the displacement current in the expanding universe, and 
we have introduced the current Jq = (Q;em/27r)8B. 

It is useful to observe that in the early universe a 3> 
H. Taking into account the expression for the Hubble 
parameter H ~ l.^Ggl^^T^ /mpi, where is the total 
number of effectively massless degrees of freedom and 
mpi is the Planck mass, we get H/a ~ 10~"'^^GeV/r 
(here we have assumed B < B^, and we have taken k ~ 1 
and 5, 10^). Note that H/a is even smaller for B > Be- 

Observing that UdI/Ij | ~ H/a, in Eq. ^ we can ne- 
glect the displacement current with respect to the Ohmic 
current. Taking the curl of Eq. Q we then get 



B 



-2HB 



V^B _^ Qfcm eV X B 

aa^ 2tt era 



(6) 



The first term in the right hand side of the above equa- 
tion describes the adiabatic dilution of the magnetic field 
due to the expansion of the universe, the second one takes 
into account the Ohmic dissipations, while the third one 
violates P and CP— symmetries, and then is responsible 
for the generation of a helical component of the magnetic 
field. Equation © describes the well known a^— dynamo 
effect (in the expanding universe) where the dynamo 
coefficient is adyn = (Q^Gm/27r)(0/(T). If the dynamo is 
"efficient" (i.e. if the the dynamo term in the right hand 
side of Eq. (0) dominates the other two) then the mag- 
netic field is exponentially amplified, while if the dynamo 
is not efficient (and neglecting any dissipative effect) the 
magnetic field is frozen into the plasma. We will see that 
this is indeed our case. 

It is useful to work in Fourier space and define the 
Fourier transform B(k,i) of the magnetic field B(x, t) 
according to B(k, i) — J (Px e~'^^'^'B{x,t), where x 
and k are co-moving coordinates and wavcnumbers, re- 
spectively. Introducing the orthonormal helicity basis 
{e+,e_,e3}, where e± = (ei ±ie2)/V2, 63 = k/fc 
(with k = |k|), and {61,62,63} form a right-handed or- 
thonormal basis, the magnetic field can be decomposed 
as B = i?+e_|_ -I- i3_e_, where i?+ and i?_ represent the 
positive and negative helicity components of B, respec- 
tively. In this basis Eq. © becomes 



B± = -2HB± 



k'^B± ^ Qfcm QkB± 



(7) 



aa" 27r aa 
It is clear from the above equation that the two helicity 
states of the magnetic field evolve differently due to the 
P and CP— violating dynamo term. The solutions of 
Eq. |7J are easily found 

S±(fc,<) -P±(fc,t,) (^)%xp(-fc2£2 ^ He), (8) 

where ti is an arbitrary time, ai — a{ti), and we have 
defined the co-moving dissipation and dynamo lengths 

il{t)= £e(0 = -^ A^-. (9) 

Jf. aa ZTT Jf. aa 



It is useful to introduce the spectra of the magnetic en- 
ergy and the magnetic helicity 

£B{k,t)=f^]\\B+\' + \B^\'), (10) 



7^B(fc,t) = a^^(|P+p-|P-'2^ 



(11) 



so that the magnetic energy and helicity are given by 

EB(t)^]^j dPxB' = j dk8B{k,t), (12) 

Hsit) = a^Jd^xA-B = JdknB{k,t). (13) 

It is clear from Eqs. 1)10(1 and pi(l that any magnetic 
field configuration satisfies the realizability condition 
\nBik,t)\ < 2a^k-^£Bik,t). Now, defining Hg'^'^(fc, t) = 
2a'^k~^£B{k,t), and inserting Eq. © in Eqs. (^01 and 
((TT|l . we get 

£B{k,t) £B(/c,ti)( — y exp(-2fc2^^)cosh(2Me), (14) 



HB{k, t) = -7^g^''(fc, U) exp(-2fc2^2) sinh(2H0), (15) 

where we have supposed that the initial helicity is null. 
We shall assume that the initial magnetic energy spec- 
trum can be represented by the following simple function 



fB(fc,i.) = ABfcPexp(-2fc24), 



(16) 



where \b and £b are constants. For k ^ l/^s, the mag- 
netic energy spectrum possesses a power law behavior, ^ 
while for large k the spectrum is suppressed exponen- 
tially in order to have a finite energy. Also, we introduce 
the so-called correlation length ^^(i) defined by 



^B{t)^a{t) 



Jdkk-^£B{k,t) 
Jdk£B{k,t) 



(17) 



The constant £b is then related to the initial correlation 
length ^B{ti) by 



r((l+p)/2) ^b{U) 
V2r{p/2) ai 



(18) 



where T{x) is the Euler gamma function. Then, £b is 
proportional to the co-moving initial correlation length. 

We now return to the evolution equation for the O- 
angle. Assuming that the back-reaction [that is the 
right-hand side of Eq. of the electromagnetic field 



^ Analyticity of the magnetic field correlator {Bi{ri)Bj {r2)) de- 
fined on a compact support for ces the spectral index p to be even 
and equal or larger than 4 l25l '). 



4 



is negligible (this will be justified a posteriori), the 
evolution of the cosmic axion field follows the stan- 
dard description Here, we just remember the 
expression for the temperature-dependent axion mass, 
ma(T) ~ 0.1m,(0){T/AQCD)-^■^ valid for ttT/A^cu > 
1, where the axion mass at zero temperature is ma(0) ~ 
0.6 (10^GeV//a) eV. The temperature at which the axion 
field starts to oscillate is usually indicated by Ti, and it 
is found by solving ma{Ti) — 3H{Ti). Its value is given 
by Ti/GeV ~ O.Qf^-^ '^A^f, where /12 = ^/(lOi^GeV) 
and A200 = AQCD/(200MeV). 

Introducing the normalized time t = t/ti, the equation 
of motion for the 0-angle becomes 



e" 



2t 



■ e' 



T^'^Q = 0, 



(19) 



where a prime indicates a derivative with respect to r. In 
Fig. 1 , we plot the solution of the above equation, where 
e,; = Q{t,) and QiU) = 0. 




FIG. 1: e(t) as a function oit/ti. 

Inserting Eqs. fTI|l and lfT5|l in Eqs. (fT^ and (fU^ we 
find, respectively 



^2 



f2 , /'2 



(l+P)/2 



X iFi 



1+P 1. ^1/2 



(20) 



HBit) = -2(27r)3a4B2 



f2 , /72 



(l+P)/2 



X iFi 



1+p 3^ £l/2 



(21) 



From Eq. H20(l we see that the dynamo is efficient (and 
then there is exponential growth of the magnetic field) if 
the dynamo length £q is greater than both the dissipa- 
tion length £d and the initial correlation length £b- 
Magnetic fields whose correlation lengths are smaller 
than the dissipation length and for which the dynamo 
is not efficient are quickly dissipated. This is easily seen 
from the evolution equation for the magnetic field, where 
the suppression factor due to dissipation is equal to the 
third term on the right-hand side of Eq. (|20(l . Because 
we are supposing that the magnetic field interacting with 
the axion is ultimately the field that we observe today in 
the galaxies and galaxy clusters, we assume that the mag- 
netic field is in the inertial regime, £d £b, during all 
its evolution. 

Since, in general, the conductivity depends on the 
value of B, to know the behavior of the dynamo length 
£q we have to solve Eq. (|20|l which is a integral equation 
for B. To this end, we re- write £q as 



fe(r) 



2tt aioi 



fir), 



(23) 



where the subscript "1" indicates that the quantity is cal- 
culated at t = ti, and we have introduced the function 
/(r) ^-JJdTT-^/^Q'/e,)iai/a), with r, = U/h. We 
have solved numerically Eq. I|20|l for f{t) taking, for pur- 
pose of simplicity, a = T(l + B/B^), which smoothly 
interpolates between the two limiting cases B <^ B^ 
and B ^ Be- In Fig. 2, we plot f{t) as a function 
of t/ti, where 6 is the solution of Eq. l(T^ . for three 
different cases of initial strength of the magnetic field. 
Here, we have parameterized the initial magnetic field as 
B, = hTf, so that we have Bi ~ 1.5 x W%{Ti/GeV)'^G. 
[If we force the magnetic energy density to be less than 
the energy density of the universe, given in the radiation 
era by p = (7rV30)g, T^, we get b < O.lg^T,)^/^ .] 




where Bi — B(ti), ifi(a, b; x) is the Kummer confluent 
hypergeometric function |26| , and we have introduced the 
root-mean-square value of the magnetic field 



B' 



'^''B(k).B(k') 



(27r)3 (27r)3 



(27r)3 



(22) 



FIG. 2: f{t) as a function of t/ti for Ti = lOGeV, Qi = 1, 
Bi = bTi, Cs(Ti) = W~''dH{T,), and p = 4. Solid line: 
b — 10^^, corresponding to Bi ~ lO^^G; Dashed line: b = 
10"'', corresponding to Bi ~ lO^'^G; Dotted line: b = 10"^", 
corresponding to Bi ~ 10^^ G. 
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We checked that, varying the initial conditions for the 
magnetic field [i.e. Ti, b, £,B{Ti), andp], the characteris- 
tic behavior of /(<) in Fig. 2, that is f{t) ~ for t < ti, 
and f{t) ~ 1 for t > ti, does not change. This is easily 
understood if we approximate (see Fig. 1) the 0-angle as 
&{t)/Qi ~ l — 9{t—ti), where 9{t) is the step-function. In 
this case, it turns out that / is a step-function centered 
at t = h, f{t) ~ e(t - ti). 

The dynamo length is small compared to the dis- 
sipative length (and then is smaller than is which we 
assumed much greater than (.d)- In fact, for t > ti and 
taking / ~ 1, we get ie/^d ^ 10~^9j/t^/^ (where we 
have taken B < Be, and then the value of io/ld is even 
smaller for magnetic fields for which B > Be). Then, 
from Eqs. l|2()|l and l|21|l . we get that in the inertial regime 
the magnetic field evolves as S oc (i.e. the field is 
frozen into the plasma) while the magnetic helicity is 



of the axion condensate evaporation. Indeed, Eq. H27|) 
can be written as 



HBit)~^2{2nYa^Btie 



(24) 



Finally, taking £q ~ (acm/27r)(0i/criai) for t > ti, we 
get that the magnetic helicity remains constant after its 
generation (say at i ~ ti) and is equal to 

Hb -b^Qi (10""G)2 kpc (25) 

for Bi < Be, corresponding to b < 10-6(GeV/Ti)2, and 



Hp 



(10-'°G)2kpc (26) 



in the case Bi > B^, or b > 10-^{GeV /Ti)^ . 
Since we expect that Qi is of order unity, and because 
Ti ~ 1 GeV, we conclude that for very strong magnetic 
fields (corresponding to 5 ~ 1), the generated magnetic 
helicity at t - ii is about \Hb \ ^ {10~^°Gy kpc. 

We can now justify a posteriori the assumption that 
the back-reaction of the electromagnetic field on the evo- 
lution of Q is negligible. ^ Integrating Eq. Q with re- 
spect to (Px, observing that Hb — ~2a'^J (PxE ■ B, and 
taking into account Eqs. Q and (|24|l . we obtain 



e + {3H + e)& + mle = 0, 



(27) 



where we have defined e = 2'na1^[ai/aY{Bl/ f^a). 
Writing Bi = bT^, and assuming fa ^ 10^^ GeV and 
B <_Bc, we get e/H ~ (e is even smaller 

for B > Be). The smallness of e justifies our previous 
assumption of neglecting the back-reaction of the elec- 
tromagnetic field on the cosmic evolution of the axion. 
We finally note that the parameter e measures the rate 



^ Below the QCD phase transition, when the axion mass is con- 
stant and equal to its zero temperature value, the fact that a 
primordial magnetic field does not interfere with the cosmic evo- 
lution of the axion has been shown in Ref. |21| . 



rf(ln Ug) 

di 



(28) 



where Ua is the number of axions per co-moving volume. 
Thus, we see that the presence of an external magnetic 
field opens a new channel for the dissipation of the ax- 
ion condensate, though this dissipation is very small (for 
other dissipation mechanisms of the axion condensate see 
|27j and references therein). 



III. CONCLUSIONS 

In this paper we have studied in detail the axion- 
magnetic field system during the period of primordial 
axion oscillations (T ^ 1 GeV). Here we summarize our 
main results and give some perspectives. 

The presence of a primordial magnetic field at the time 
when the cosmological axion oscillations begin amplifies 
the axion decay probability and, contemporaneously, cat- 
alyzes the photon production (enlarging the magnetic en- 
ergy itself). The effect of this process on the axion field 
is minimal, since the number of decaying axions is rather 
small Aua/ua ~ 10~^°. A simple consequence is that the 
cosmological bound on the axion mass is preserved, even 
in the presence of a very strong external magnetic field. 

However, the axion oscillations leave a signature on 
the magnetic field itself, in the form of magnetic helic- 
ity. This is due to the different probability for creat- 
ing photons of different helicities. Thus the axion os- 
cillations, interacting with a magnetic field, originate a 
macroscopic CP— odd state. The amount of helicity pro- 
duced depends on the initial intensity of the magnetic 
field. At late times, the turbulence of the primordial 
plasma should be taken into account in considering the 
evolution of both magnetic energy and helicity densi- 
ties. However, it is now believed that helicity is a quasi- 
conserved quantity in magnetohydrodynamic turbulence 
[^H^. Hence, the helicity generated in the primordial 
universe should survive until today without changing its 
major properties, and could be important for the dy- 
namo mechanism of magnetic fields operating in galaxies 
and clusters of galaxies (see, e.g., Ref. 0). Moreover, 
the presence of a helical component may speed up the 
growth of the magnetic field correlation length during its 
evolution in the primordial universe through the so-called 
inverse cascade of magnetic energy l29j , leading today to 
magnetic fields on larger scales [a ISfl ■ 

In principle, a helical magnetic field could leave pe- 
culiar imprints on the Cosmic Microwave Background 
Radiation (CMBR). Unfortunately, the maximal helicity 
producible in our mechanism, \Hb\ — (10~^°G)^ kpc, is 
much smaller than that detectible in near future CMBR 
experiments, which is of order of (10~^G)^Mpc [sj. 
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However, from the above discussion it emerges a strong 
relation between axion and magnetic helicity. We believe 
that this subject deserves further investigations. 
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